INFLATING A CHAIN OF X-RAY DEFICIENT BUBBLES BY A SINGLE 

JET ACTIVITY EPISODE 



Michael Refaelovich 1 and Noam Soker 1 



ABSTRACT 

We show that a continuous jet with time-independent launching properties can inflate a 
chain of close and overlapping X-ray deficient bubbles. Using the numerical code PLUTO we run 
2.5D (i.e. spherical coordinate system with cylindrical symmetry) hydrodynamic simulations and 
study the interaction of the jets with the intra-cluster medium (ICM). A key process is vortex 
fragmentation due to several mechanisms, including vortex-shedding and Kelvin-Helmholtz (KH) 
instabilities. Our results can account for the structure of two opposite chains of close bubbles 
as observed in the galaxy cluster Hydra A. Our results imply that the presence of multiple pairs 
of bubbles does not necessarily imply several jet-launching episodes. This finding might have 
implications to feedback mechanisms operating by jets. 



1. INTRODUCTION 



Bubbles (cavities) devoid of X-ray emission, mostly as opposite pairs, are observed in a la rge fraction 
of co oling flow clusters and groups o f galaxies, as well as in cooli ng flow elliptical galaxies (e.g.jDong et al 



01 co oling now clusters and groups 01 galaxies, as well as m cooli ng now elliptical galaxies (e.g.JJJong et al. 
l20ld ). Examples inclu de Hydra A ([Wise et al.ll2007|), RBS7 97 ([Doria et al.ll2012l: ICavagnolo et al.ll201lh . 



Perseus JFabian et al.l l201lL A bell 2052 falanton et a lJ boilh. NGC 6338 (IPandge et alJIgQl^. NGC 5044 
(jPavid et al.ll201ll), NGC 4636 feakhet alfcoOQh . NGC 5044 (|Pavid et allboooh . HCG 62 (|Gitti et alboioh 
and NGC 5044 ([Gastaldello et aL 2009[ ). These bubbles are inflated by jets launched from the central active 
galactic nuclei (AGN), as evident by the radio emission that f i lls most bubbles. Th e se jets and cavities 
heat the i ntracluster medium (ICM; e.g.. lO'Sullivan et al.l 12011c IGaspari et al. |[2012a,b ; IBirzan et al.l 1201 ll 



Gitti et al. I l2012l for recent papers and references therein), and maintain a negative feedback mechanism with 



the cooling gas ( Binney fc Taborl 1995 : Far age et al. 2012 ). 



Wide bubbles very close to the origin of the jets (the AGN), e.g., as in Abell 2052, that are termed 'fat 
bubbles', can be inflated by jets that do not penetrate through the ICM. Instead, they deposit their energy 
relatively close to their origin and inflate the fat bubbles. Slow massive wide (SMW) jets can inflate the 
fat bubbles that are ob served in many cooling flows, in clusters, groups of galaxies, and in elliptical galaxies 



([Sternberg et al. II2QQ7I ). The same basic physics that prevents wide jets from pe netrating through the ICM 
holds for precessing jets ([Sternberg fc Soker 2008a; Falceta-Goncalves et al. 2010) or a relative mo tion of the 



jets to the medium ( Briiggen et al. 20071 : ISoker 2009 : 



Morsonv et al. 



2010l : lMendvgral et al.ll2012l ). If the jet 



penetrate to a too large distance, then no bubble s are formed, while in intermediate cases elongated and /or 



detached from the center bubbles a re formed (e.g., Basson fc Alexander! 2 003 [ 



Mendvgral et al.ll201lD . 



2004; 



Heinz et al. 



_ Omnia et al 

2006 : IVernaleo fc Reynold] I2OO6I : lAlouani Bibi et al.l l2007i Isternberg et al 1 120071 : lO'Neill fc Jones! [201 
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Vortices inside the bubbles and in their surro undings play major roles in the formation of bubbles, 
their evolution, and their interaction with the ICM. I Omnia et al. ( 2004h find that a turbulent vortex trails 
each cavity, and tha t this vortex contains a significant quantity of entrained and uplifted material (also 
Roediger et al J 120071). Jet-excited shocks that interact with older bubbles can excite vortices that dissipate 
energy to the ICM ( Friedman et al.|[2012 ). The vortices cause semi-periodic change s in the bubble properties , 
such as its boundary. This causes a single bubble to excite several sound waves ([Sternberg fc Sokerll200^ . 
Vortices play a majo r role in mixing the shocked jets' material with the ICM, hence heating the ICM gas 
(IGilkis fc Sokerll2012J). The role of vor tices in the interaction of the jet with the ICM was studied before ( 
e.g., iNormanl 19961 : iMizuta et al. 2004 ). However, no special emphasize on bubbles formation was done. 



In some ca ses two opposite c hains of bubbles that are close to each other, and even overlap, are observe d 
as in Hydra A ([Wise et al.l 120071 ). and in two bubbles in the galaxy group NGC 5813 ([Randall et al.l 120111 ). 
These chains were usually attributed to several episodes of jet activity. Using a 2.5D hydrodynamical code 
(section [2} we study the formation and evolution of chain of bubbles from one continuous jet activity episode 
(section [3]). Our short summary is in section HJ 



2. NUMERICAL SET UP 



We use the multidimensional hydrodynamic PLUTO code ([Mignone et al.ll2007l ) for all our simulations. 
The simulations are 2.5D in the sense that we use a spherical coordinate system, but impose cylindrical 
symmetry, hence calculating the flow with a 2D polar grid (r, 6). There is also a mirror symmetry on the 
plane = tt/2. There are 256 divisions in the azimuthal direction (0 < < ^) and 257 divisions in the radial 
direction, geometrically stretched in the range (1 < r < 400 kpc). Our flow is non-relativistic. Gravity is 
included, but as the simulations last for a time shorter than the radiative cooling time we neglect radiative 
cooling. The exact scheme that we use includes third-order Runge-Kutta time-stepping integration (two 
predictors), parabolic interpolation and no dimensional splitting. 

The boundary conditions are reflective on both 'angular' boundaries: the symmetry axis 6 = and the 
equatorial plane 6 = 7r/2, and an outflow condition on the far radius (outer radial boundary). The jet is 
injected by using the boundary conditions at r = 1 kpc. The jet is injected within an angle < < a 
- where a is the half-opening angle. At other angles on the inner radius we impose reflective boundary 
conditions. 



The initial density profile is based on that of Hydra A as reported bv iMcNamara et aL ( 2000 ). In the 



inner 10 kpc we approximate the profile by a constant density of n e = 0.06 cm -3 , where n e is the electron 
density, while at r > 10 kpc the density drops as n e (r) oc r _1 . The initial temperature of the gas was 
constant across the entire domain To = 3.5 keV. The gravity is calculated from the initial pressure and 
density profiles by assuming hydrostatic equilibrium g = — • 

We performed many simulations spanning a large volume of the parameters space, e.g., various opening 
angles, jet powers, and jet velocities. The standard case is taken to have the following parameters: Total 
two jets power of P^j = 4 x 10 45 erg s _1 , a jet Mach number of 15 relative to the ambient medium, a jet 
velocity of Vj = 1.33 x 10 4 km s _1 , hence the mass outflow from the two opposite jets (here we calculate 
only one jet) is M^j — 7OM yr _1 , and a half-opening angle of a = 30°. 
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3. RESULTS 



3.1. Morphology 

In figured] we present our standard run at time t = 100 Myr. We focus our study on the vortices that 
can be clearly seen along the shocked jet's material, i.e., the cocoon. Between the large vortices we can see 
denser ICM gas segments flowing toward the axis. In our simulations the axi- symmetry implies that each 
such segment is actually a torus. These segments are marked by "KH instability". There is a bow shock 
running into the ICM, the forward shock. Most of the kinetic energy of the jet is dissipated in the reverse 
shock, as marked on the figure. There is a contact discontinuity surface between the shocked ICM and the 
shocked jet material. It can be located where the large temperature and density gradients are. 

The primary source of vorticity in the simulation is the back-flow of the shocked jet material reflected 
from the shock front. As may be seen in figure [2j a large vortex immediately behind the jet's head, the 
"primary vortex" , exists at all times, and it creates a cascade of smaller vortices. 



Morphologically, these vortices can be identified with bubbles (see also ISternberg fc Sokerll2008bf ). Sev- 
eral process, as we now detail, influence the evolution of these vortices, and hence the bubbles' morphology. 



• Kelvin- Helmholtz (KH) instability. The back- flow of the lower density cocoon material inside the 
vortex relative to the denser shocked ICM material forms a shear flow that is prone to the KH instability. 
The KH instability creates the wavy morphology, namely, it fragments large bubbles to smaller ones. 

• Rayleigh- Taylor (RT) instability. The denser shocked ICM is on top of the much l ower density gas of 
the b ubbles. This RT instability does not exist in the front of the outer bubble ([Sternberg fc Soker 



2008bf ). but it might amplify the fragmentation formed by the KH instability in the tail. 



• Vortex S hedding. Vortex shedding occurs mainly on trailing vortices and forms disconnected vortices 
(|Normanlll996[ ). 

• Jet flapping. The jet is collimated by the pressure of the earlier shocked jet material and the shocked 
ICM. Due to the instabilities listed above and the stochastic flow the collimating pressure is not 
constant. As a result of that the jet's width changes. In 3D flow the jet will jitter around the 
symmetry axis. This motion is termed flapping. This flapping occurs in semi-periodic manner on a 
typical time scale of jet-width sound crossing time. This flapping has two effects. First, the jitter of 
the jet slows down the 'drilling' of the jet through the ICM. Second, the flapping seeds semi-periodic 
perturbation on the boundary of the primary bubble which in turn amplified by the KH instability 
and at least at early stages, creates vortex-shedding. The flapping is the subject of a future paper. 



The type of flow simula ted here is known to be KH- unstable, with formation of m ore than one vortex in 
the non-linear regime (e.g., Norman 1996 : Mizuta et al. 2004 : Mathews fc Guol 2012 ). Our new emphasize 



is on the role of these large vortices in forming and determining the properties of observed X-ray deficient 
bubbles (cavities). Although magnetic fields might change the exact morphology of the vortices formed by 
the KH instability, they cannot suppress them compl etely. If the magnetic fields are entangled on small 
scales, as is expected in many cases (e.g., ISokerl 120101 ). then the magnetic field can locally suppress small 
scale KH-instability modes. However, they cannot suppress the large scale modes that are considered here. 

In the opposite limit to small scale fields, the ICM magnetic field lines are stretch ed by the jet and 
exert tension on two opposite sides of the shocked jet's material (see, e.g., figure 3 in LSokerll2010l ). If the 
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Fig. 1. — Results of the standard run at t = 100 Myr. Axes are in kpc. The horizontal boundary (labelled 
Z) is the symmetry axis of the two jets. The left vertical boundary of the figure (labelled R) is taken 
along the mirror-symmetry plane of the flow z = 0; it is termed the equatorial plane. The total two jets 



power is P2j = 4 x 10 erg s , the initial jet velocity is Vj 



- 1.33 x 10 4 km s _1 , the mass outflow from 
the two opposite jets (here we calculate only one jet) is M^j = 7OM yr _1 , and the jet half-opening angle 
is a = 30°. The jet is injected at r = 1 kpc, a region that is not well resolved in the figure. We mark 
several features that are discussed in the text. |(a)| The temperature map of the standard run in logarithmic 
scale and units of K. |(b)|The electron density in logarithmic scale and units of cm -3 , (c) Contours of the 



electron density in logarithmic scale and velocity arrows. The velocity arrows are divided into 4 velocity 
bins: 500 - 1000 km s" 1 , 1000 - 5000 km s~\ 5000 - 10 4 km s _1 and above 10 4 km s -1 - from shortest to 
longest. Velocities below 500 km s _1 do not appear in the figure. The fragmented vortices, that are at the 
heart of the present study, are marked. Vortex shedding region is also marked. 
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magnetic fields are strong enough they can suppress even the large scale KH-instability modes, but only along 
the field lines. Namely, only in two opposite quarters around the jet. The two other quarters are prone to 
the KH instability and bubbles will be formed. I n any case, strong magnetic field that can suppress large 
scale KH-instability modes are not easy to obtain. llVIathews fc GuJ ( 2012 ). for example, did not succeed in 
suppressing the KH instability in their MHD simulations. 



3.2. Chains of cavities 



The aforementioned mechanisms of vortex formation and fragmentation lead to the formation of a chain 
of vortices along the jet axis (figure |l(c)| ). To compare with observation we integrated the quantity along 



the line of sight perpendicular to the symmetry axis. To a good approximation n\ is proportional to the 
observed X-ray emission. To emphasize the cavities we a pplied Gaussian u nsharp filter to the integrated 



map, with filter parameters of cr coarse = 30 and <jfi ne = 10 (|Dong et al.ll2010f ). The results are presented in 
figure 03 We can identify two large bubbles and a few smaller ones closer to the center at each side of the 
equatorial plane. The center of the large bubbles are located approximately at coordinates (r, z) = (0, ±150) 
and (0, ±210). The denser medium between the large bubbles which creates a brighter X-ray emission and 
separates the bubbles is seen at z = ±130 and z = ±175. 

Our results can account for the structure of two op posite chains of bu bbles that are close to each other, 



and even overlap, as observed for example in Hydra A (|Wise et al.l 120071 ) and in two bubbles of NGC 5813 



([Randall et al.ll201lf ). The inner smaller bubbles that are seen in our simulations are likely to be smeared in 
observations due to departure from pure axi-symmetry, inclination, and low resolution, and appear as one 
bubble. However, under different numerical parameters, more than two large bubbles on one side can be 
obtained. 

We now emphasize two important points. First, we note that because of the highly turbulent nature 
of the flow we expect the simulations not to be convergent in terms of fine structural details, i.e. small 
vortices and thin filaments. Moreover, small differences in the initial parameters or numerical noise may 
result in macroscopic differences due to these instabilities. Such differences may include the exact position of 
the forward shock, the number of macroscopic vortices and whether the jet penetrates through the primary 
vortex. The total volume of the bubbles is determined by energy considerations, and does not depend on the 
these instabilities. As well, the general large-scale morphology does not depend on the instabilities. Indeed, 
we find that chains of bubbles are formed for a non- negligible set of jet parameters (see section l3T3j) . 

Second, in order to appear like a chain of bubbles in the surface brightness map, the vortices shouldn't 
necessarily be completely disconnected. It is sufficient that a relatively small ICM segment starts to fragment 
a bubble into two regions, for the structure to appear like two adjacent bubbles in the X-ray image. In fact, 
in most our simulations the bubbles are at least partially connected by the jet. Such dense ICM segments 
are regularly formed by vortices in our simulations, and hence a chain of bubbles is not a rare outcome. 



3.3. Parameter space survey 

We examine the influence of the jet's half opening angle a and of its velocity Vj on the appearance of 
the chain of bubbles viewed perpendicular to the jets' axis. We have studied the range 15° < a < 60° and 
7 x 10 3 km s _1 < vj < 3 x 10 4 km s _1 . In figure [4] we present the calculated X-ray map of a subset of the 
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runs after applying the unsharp filter (as explained in section |3T2|) . 

In most of these synthetic images bright X-ray thin segments divide the X-ray deficient volume into 
several apparent bubbles. These thin segments are the KH-filaments of the ICM pushing between the 
bubbles (the area marked by 'KH-instability' in figure |l(a)| ) • In reality, inclination and departure from pure 
axisymmetry will make the segments less prominent. Therefore, in observations only cases with high segment 
to bubble contrast will appear like a chain of bubbles. ^From figure 2] we find that these cases are typical 
for mid-range opening angles (a = 30°, figures |4(d)|4(e)|4(f)1 ), but with some dependence on the jet velocity. 



In cases with slower jets (vj < 1.3 x 10 4 km s _1 ) and larger opening angles (a > 45°) as in figures [4(g)] and 
|4(h)| they will mostly appear like a single bubble. For very fast and narrow jets, the jet escapes without 
creating bubbles (figures |4 (a) |4(b)|4(c) ). 



4. SUMMARY 



We showed that a combination of various vortex fragmentation mechanisms, mainly KH-instability and 
vortex-shedding, can explain chains of close and over lapping X-ray cavitie s (bubbles) observed in clusters of 
galaxies. We used the hydrodynamic PLUTO code (jMignone et al.l 120071 ) in 2.5D, i.e. spherical coordinate 
system with cylindrical symmetry. The general flow structure is presented in figure [TJ where two large 
bubbles can clearly be identified. A primary (front) vortex exists during the entire simulation with one or 
two large trailing vortices and smaller vortices shedded further downstream (figure [2} . The large vortices are 
interpreted as the X-ray deficient bubbles seen in observations, as demonstrated in figure [3l Although chain 
of bubbles are formed in a wide range of parameters, for some parameters, like for wide jets, the bright X-ray 
segments separating bubbles are too faint and only one bubble is formed for each jet launching episode. This 
is demonstrated in figure HJ 

We find the formation of many vortices by the KH-instability and the vortex shedding process to be 
very interesting and significant in determining the structure and evolution of bubbles. It may also play an 
importan t role by dredging m aterial up from the inner regions outward and by dissipating heat to the inner 



region (jGilkis fc Sokerll2012h . In a future study we will explore t he entrainment and lifti ng of ICM by the 
chain of vortices, a process that is observed to occur in Hydra A (jKirkpatrick et al.l l2009). 



Basically we have shown that one jet activity e pisode can form two opposi te chains of close bu bbles as 
observed for example in Hydra A ([Wise et al.ll2007[ ) and possibly in NGC 5813 ([Randall et al.ll2Qllh . Large 
separations between bubbles, on the other hand, probably require multiple jet activity episodes. 

We thank an anonymous referee for helpful comments. This research was supported by the Asher Fund 
for Space Research at the Technion, and the Israel Science foundation. 
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Fig. 2. — Density maps of the standard run at four times. The color scaling is of logn e (cm -3 ), where n e 
is the electron density. A primary vortex right behind the jet's head and another vortex trailing it exist 
during the evolution. At the t = 10 Myr picture the flapping of the jet may be seen as a perturbation on 
the separating layer between the expanding jet and the cocoon. 
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Fig. 3. — Synthetic X-ray image of the standard simulation viewed from the equatorial plane direction 
(inclination of 90°) at t = 100 Myr, performed by integrating v? e along the line of sight, (a) Original X-ray 
image. |(b)| Unsharp filter applied to the image at |(a)| Thin lines are the 1 contours of the quotient image. 
The noisy contours which are not on the major axis are due to sound waves. 
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(a) a = 15°, ^ = 6.5 x 10 3 km s" 1 (b) a = 15°, vj = 1.3 x 10 4 km s" 1 (c) a = 15°, ^ = 2.7 x 10 4 km s" 1 
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(g) a = 45°, vj = 7 x 10 3 km s -1 (h) a = 45°, vj = 1.3 x 10 4 km s -1 (i) a = 45°, vj = 2.7 x 10 4 km s -1 

Fig. 4. — Synthetic X-ray images as in figure [3j but for different values of the jet's half opening angle a 
and velocity Vj, as specified in each panel. Panel [(e)] is our standard case. The rest of the parameters are 
identical to those of the standard run (see section [2]). 



